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I. Introduction 


Some experimental spacecraft use superconducting sensors for gyro read-out 
and so must be maintained at a very low temperature. The boil-off from the 
cryogenic liquid used to cool the sensors can also be used, as the Gravity Probe- 
B (GP-B) spacecraft does, as propellant to maintain attitude control and 
drag-free operation of the spacecraft 1 . The cryogenic liquid for such spacecraft 
Is, however, susceptible to both slosh- like motion and non- axisymme trie 
configurations under the influence of various kinds of gravity jitter and gravity 
gradient accelerations. Hence, it is important to quantify the magnitude of the 
liquid- induced perturbations on the spacecraft. Dynamic behavior of liquids in 
a moving container under reduced gravity environment have been widely discussed 
by Abramson 2 . Here we use the example of the GP-B to investigate such 
perturbations by numerical simulations. For this spacecraft disturbances can be 
imposed on the liquid by atmospheric drag, spacecraft attitude control maneuvers, 
and the Earth's gravity gradient. More generally, onboard machinery vibrations 
and crew motion can also create disturbances. Recent studies 3-7 suggest that 
high frequency disturbances are relatively unimportant in causing liquid motions 
in comparison to low frequency ones. The results presented here confirm this 
conclusion . 

After an initial calibration period, the GP-B spacecraft rotates in orbit 
at 0.1 rpm about the tank symmetry axis. For this rotation rate, the equilibrium 
liquid free surface shape is a "doughnut" configuration for all residual gravity 
levels of 10" 6 g o or less, as shown by experiments 8 and by numerical simulations 4-7 ; 
furthermore, the superfluid behavior of the 1.8 K liquid helium used in GP-B 
eliminates temperature gradients and therefore such effects as Marangoni 
convection do not have to be considered. Classical fluid dynamics theory is used 
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as the basis of the numerical simulations here, since Mason's experiments 9 show 
that the theory is applicable for cryogenic liquid helium in large containers 10 . 
To study liquid responses to various disturbances , we investigate and simulate 
three levels of gravity jitter (1CT 6 , 10" 7 , and 10~ 8 g o ) each at three predominant 
frequencies (0.1, 1.0, and 10 Hz), combined with a gravity gradient appropriate 
for the GP-B orbit 11 " 13 . In this study, dynamical evolution of sloshing dynamics 
excited fluid forces and torque fluctuations exerted on the dewar container 
driven by the combined gravity gradient and jitter accelerations are also 
investigated and simulated. 

II. Mathematical Model of the Simulation of Sloshing 
Dynamics Modulated Bubble Oscillations 
Consider a closed circular dewar of an inner radius a x and an outer radius 
a 2 ; top and bottom dome of radius b; with height L, which is partially filled 
with a cryogenic liquid helium and the rest of the ullage is filled with a helium 
vapor. Density and viscosity of liquid helium and helium vapor are p L , /i L , p v , 
and /Xy, respectively. Let us use cylindrical coordinates (r, 0 , z) , with 

corresponding velocity components (u, v, w) , corresponding components of gravity 
gradient acceleration (a gg r , a gg B , a gg z ) , and corresponding jitter acceleration 
( a sj,r» a gj ,$* a gj,z)‘ Th e whole fluid system is spinning in the z-axis with an 
angular velocity of and an angular acceleration of a. The governing equations 
for non- inertia frame spacecraft bound coordinates are illustrated in our recent 
studies 14,15 , and will not be repeated in this paper. In these formulations, 2pc*)u 
and 2 pwv denote the Coriolis forces; pru> 2 , the centrifugal force; and pra, the 
angular acceleration of spinning motion. 

In the computation of fluid forces, moment, viscous stress and angular 
momentum acting on the container wall of the spacecraft, one has to consider 
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those forces and moment in the inertia frame rather than the non- inertia frame. 
For the case of spinning motion in the z-axis, one has to transform those vectors 
from the non- inertia frame to the inertia frame as follows: 
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where the prime symbol denotes vectors in the inertia frame while those 
parameters without the prime symbol indicate vectors in the non- inertia frame. 

In order to solve sloshing dynamic problems of liquid propellent systems 
In orbital spacecraft under microgravity environment, one has to solve the 
governing non- inertia frame equations 14 - 15 accompanied by a set of initial and 
boundary conditions. A detailed illustration of these initial and boundary 
conditions concerning sloshing dynamics of fluid systems in microgravity were 
precisely given in our recent studies 15 - 16 and will not be repeated in this paper. 
In this study, to show a realistic example of the simulation, a full scale GP-B 
spacecraft propellant dewar tank with an outer radius of 78 cm and an inner 
radius of 13.8 cm, top and bottom radius of 110 cm and a height of 162 cm has 
been used in the numerical simulation. The dewar tank is 80% filled with 
cryogenic liquid helium and the rest of the ullage is filled with helium vapor. 
The temperature of cryogenic helium is 1.8 K. In this study the following data 
was used: liquid helium density - 0.146 g/cm 3 , helium vapor density - 0.00147 
g/cm 3 , fluid pressure - 1.66 x 10 4 dyne/cm 2 , surface tension coefficient at the 
interface between liquid helium and helium vapor - 0.353 dyne/cm, liquid helium 
viscosity coefficient - 9.61 x 10' 5 cm 2 /s; and contact angle - 0°. The initial 
interface profiles of the liquid-vapor-solid phase used in this study have been 
given explicitly through the steady state computations made by Hung and Leslie 17 
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and Hung et al. 4,5,17 , which were checked by the experiments carried out by 
Leslie 8 . 

A staggered grid for the velocity components is used in this computer 
program. The marker-and-cell method of studying fluid flows along a free surface 
is adopted 18-21 . The formulation for this method is valid for any arbitrary 
interface location between the grid points and is not limited to middle point 
interfaces 20 . 

There are several methods to be used to solve fluid flow simulation 
problems in a free surface. Explicit, semi -explicit and implicit schemes are the 
most commonly used methods. Mixing methods for these schemes are employed in 
this study. A near semi- implicit PCMI method (Predictor-Corrector Multiple- 
Interaction Technique) 22 is used to simulate velocity profiles, and Gauss- 
Seidel 22 , or Conjugate Gradient 2 * methods are employed to simulate pressure 
iteration, while Successive Over-Relaxation method 25 is used to simulate free 
surface configurations through iteration. It has been shown that near semi- 
implicit scheme for PCMI method is more reliable than that of the explicit 
scheme, and also more convenient and computational time saving than that of the 
implicit scheme 26 . Second order accuracy of derivatives are adopted in the 
finite-different equations while the truncation errors for the interior of time 
and spatial steps are of (At 2 ) and (Ar 2 , A 6 2 1 Az 2 ) , respectively. The convergence 
criteria for velocity, pressure and liquid volume are less than 10 -5 . Some of 
the steady state computational results are compared with the experimental 
observations carried over by Leslie 8 in free-falling aircraft (KC-135) with 
excellent agreement. 

III. Simulation of Bubble Oscillations Driven by 
Gravity Gradient and Jitter Accelerations 
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With an example of the GP-B Spacecraft, which is an Earth satellite 
orbiting at 650 km altitude directly over the poles, the orbit period, r Q can be 
computed from following expression: 


*0 = 2 * 


Rl ' 2 


R b 9o 2 


( 2 ) 


where R E denotes radius of Earth (— 6373 km); R c , radius of circular orbit (— R E 
+ h - 7023 km); h, orbit altitude (- 650 km); and g ot Earth gravity acceleration 
(" 9.81 m /s 2 ). For the case of GP-B, the orbit period r Q — 97.6 min, and orbit 
rate n - 2n /r 0 - 1.07 x 1CT 3 rad / s . 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 12 ' 27 


- -n 2 [3(f c -d) f c -d] (3) 

where a gg denotes gravity gradient acceleration vector; 3, vector from the fluid 
element to the spacecraft geometric center; r c , unit vector from the spacecraft 
geometric center to the center of the Earth; and n, the orbit rate. 

For the case of the GP-B Spacecraft, it is assumed that the gravity exerted 
on the mass center of the spacecraft orbiting around the Earth on its specified 
or bit is zero. In other words, all the gravity acceleration exerted on the 
spacecraft is nothing but the gravity gradient acceleration which is defined in 
Equation (3) . 

Figure 1 illustrates the geometrical relationship of the parameters shown 
in Equation (3). Let us consider the fluid element of interests, m, located at 
(r, 6, z) in cylindrical coordinates and at (x, y, z) in cartesian coordinates. 
As |d| is much smaller than the distance between the location of the GP-B 
Spacecraft mass center to the center of the Earth, r c through the GP-B geometric 
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center and r c through the fluid element, m, is basically the same. 

The gravity gradient acceleration acted on the fluid element located at (r, 
8 t z) , in non-inertia frame coordinates is shown as follows 14,15 : 
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3 [ ( z-Lf 2 ) cosij^-rcos (0+Q t) *sim|J E ] sinijr^+rcos (0+u> t) 

rsin(0+o>t) (4) 

-3 [ (z-L/2) cosijf^-rcos (0+co t) -siruj/^] costy s + (z-L/2) 

where ^ denotes the azimuth angle of Earth toward the location of spacecraft 

geometric center, or angle between spacecraft rotating axis and radial direction 

from spacecraft geometric center to center of Earth; n, the orbit rate; and t, 

the time of spacecraft orbiting around the Earth. 

At time t — 0, the rotating axis of spacecraft is aligned with the radial 

direction of Earth center to the spacecraft mass center. Assuming spacecraft 

axis of rotating is linearly turning around 0° to 360° in 1200 s when spacecraft 

is orbiting around the Earth. Thus, can be defined as 

^=2-|t (5) 

where r is spacecraft gravity turn-around time (- 1200 s in this case); and t is 
the time measured from the instant when the direction of spacecraft rotating axis 
is aligned with the radial direction of spacecraft mass center to the center of 
the Earth. 

Fluctuations of residual gravity due to gravity jitter acceleration is 
modeled by following equation 11 : 
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( 6 ) 


9 r =9 r B [l + -jSin(2jifC) ] 

where g B denotes the background gravity environment, and f(Hz) stands for the 
frequency of gravity- j itters . 

In this study, three ranges of gravity jitter background with gravity of 
10~ 6 , 10' 7 and 1CT 8 g Q due to different intensity of spacecraft atmospheric drag 
and rotating speed of 0.1 rpm for normal GP-B spacecraft operation, and three 
frequencies of gravity jitter with the values of 0.1, 1.0 and 10 Hz have been 
considered for the investigation of the oscillations of the liquid-vapor 
interface. Components of gravity jitter acceleration in non- inertia coordinate 
system shall be given by 11 ' 14 ** 15 

£gj= ( a gj,z> a gj.Q’ a g).z ) = tesiniJrgCos (0+0) t) , -gsimlr^in ( 0 +<o t) , -gcosijr^] 

( 7 ) 

The characteris tics of gravity gradient and j itter accelerations, shown in 
Equations (4) and (7), respectively, are quite different. Gravity gradient 
acceleration indicates that acceleration acted on any fluid mass inside the 
container increases two units of acceleration per unit of distance measured from 
the mass center of the container (point 0 & at Figure 1) to the location of the 
fluid mass parallel along the radial axis from the mass center of the container 
to the center of the Earth (parallel to unit vector r c shown in Figure 1) while 
the acceleration acted on the fluid mass decreases one unit of acceleration per 
unit of the shortest distance measured from the location of the fluid mass to the 
radial axis along the vector from the mass center of the container to the center 
of the Earth 12 . To give an example, Figure 2 shows the time variation of gravity 
gradient accelerations for a turn-around period of 1200s with container rotating 
speed of 0.1 rpm for components along (x, y, z) directions acted on the fluid 
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mass located at (r, 9 , z) - (60 cm, tt/ 4 , 40 cm). As the magnitude and direction 
of gravity gradient acceleration acted on each fluid mass is strongly dependent 
upon how far the location of the fluid mass deviates from the mass center of the 
container measured along the axis parallel to the vector r c which varies with 
respect to time, it can be anticipated that gravity gradient acceleration acted 
on the fluid mass is different for fluid mass at different locations in the 
container. Figure 2 shows that the magnitude of gravity gradient acceleration 
is on the order of 1CT 7 g^ On the contrary to gravity gradient acceleration 
which acts different acceleration on the fluid mass at different locations in the 
container, gravity jitter acceleration drives the same acceleration on the fluid 
mass at different locations in the container. Figure 3 shows the time variation 
of gravity jitter accelerations for a turn-around period of 1200 s with a 
container rotating speed of 0.1 rpm and a jitter frequency of 0.1 Hz for 
components along (x, y, z) directions acted on the fluid mass everywhere in the 
container. It is shown in Figure 2 that gravity gradient accelerations are on 
the order of 10 7 g Q for the spacecraft orbit chosen in this example, three 
ranges of background gravity, 10 6 , 10~ 7 and 10~ 8 ^ for gravity jitter 
accelerations correspond to the magnitude of acceleration higher, equal, and 
lower, respectively, than that of the gravity gradient acceleration acted on the 
fluid system of the spacecraft. 

The equilibrium shape of the liquid-vapor interface for a rotating dewar 
under a residual gravity environment below 10~ 6 g o and rotating speed of 0.1 rpm 
is very much alike and is a doughnut configuration with a near circular kidney- 
shape cross -section based on the computation of the numerical algorithm developed 
in our earlier studies' 1 . Figure 4(A) shows the initial shape of the interface 
in the r-z plane at 6 - 0° and 180°; Figure 4(B) shows the initial profile of the 
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liquid-vapor interface in the r-z plane at 9 - 90° and 270°; Figure 4(C) shows 
the initial profile of the liquid- vapor interface in the r-0 plane at height z 
- 108 cm; and Figure 4(D) shows the initial profile of three-dimensional liquid- 
vapor interface. Becuase of page limitation, the time evolution of three- 
dimensional bubble oscillations only will be illustrated in this paper. 

In this study, combined gravity gradient and three ranges of jitter (10~ 6 , 
10~ 7 , and 1CT 0 g 0 ) accelerations acting on the liquid-vapor interface oscillations 
have been simulated and investigated. 

(Ill -A) Combined Gravity Gradient and 10" 8 g 0 Jitter 
Accelerations Acting on Bubble Oscillations 
In this case, combined effects of gravity gradient and jitter accelerations 
are completely dominated by the gravity gradient acceleration when the 1CT 8 
and less gravity jitter accelerations are applied simultaneously with gravity 
gradient acceleration (which is on the order of 10" 7 g 0 as shown in Figure 2) for 
the excitation of slosh waves along liquid-vapor interface oscillations. 

Figure 5 shows the time sequence evolution of the bubble profiles driven 
by these combined accelerations. For the convenience of comparison, figures of 
liquid-vapor interface profiles with the same values of the time sequences chosen 
throughout this paper are at time t - 191, 354, 380, 431, 503, 603, 825, 980, 
995, 1050, 1080 and 1200 s. It clearly shows that there are a series of 
asymmetric oscillations excited along the surface of liquid-vapor interface 
driven by asymmetric gravity gradient- dominated acceleration. 

Careful examination for the case of gravity gradient-dominated 
acceleration- effected liquid-vapor interface oscillations , gravity gradient 
acceleration shown in Equation (3) , indicate that there are greater negative 
components of acceleration longitudinal to the direction aligned with spacecraft 
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mass center to the Earth center, and smaller positive components of acceleration 
transverse to this direction. As we indicated in Equation (5), angle ^ varies 
with time. This phenomena shows that gravity gradient acceleration exerted on 
the spacecraft is equivalent to the combination of time dependent force with 
turn-around direction and torsional moment acting on the spacecraft when it is 
orbiting around the Earth. It shows in Figures 4 and 5 that the deformation of 
the bubble is created by an asymmetric torsional moment with a twisting force. 

Figure 5 shows the time sequence evolution of the three-dimensional 
dynamical behavior of the bubble oscillations driven by gravity gradient- 
dominated acceleration. Left-side up and right-side down oscillations of bubble 
in the vertical cross-section of doughnut profiles indicate a very important 
characteristics of gravity gradient-dominated acceleration which produce a 
combination of time dependent equivalent- torsional moment and twisting force 
acting on the spacecraft when it is orbiting around the Earth. 

( II I - B) Combined Gravity Gradient and 10‘ 7 g 0 Jitter 
Accelerations Acting on Bubble Oscillations 

In this section, combined effects of accelerations are equally weighted by 
both the gravity gradient and jitter accelerations when 10~ 7 gQ gravity jitter 
acceleration is jointly applied simultaneously with gravity gradient acceleration 
(which is on the order of 10' 7 g 0 as shown in Figure 2) for the excitation of 
slosh waves along the liquid- vapor interface oscillations . Figure 6 shows the 
time sequence evolution of dynamical behavior of the bubble oscillations driven 
by these combined accelerations. Similar to Figure 5, same values of the time 
sequences are chosen. It also shows that there are a series of asymmetric 
oscillations excited along the surface of liquid-vapor interface driven by 
asymmetric equally weighted gravity gradient and jitter accelerations. 
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As we recorded in Figure 5, the effect of gravity gradient -dominated 
acceleration exerted on the spacecraft is equivalent to the combination of time- 
dependent force with turn-around direction of torsional moment acting on the 
spacecraft when it is orbiting around the Earth. On the other hand, it will be 
show later in Section (III-C) for the case of gravity jitter-dominated 
acceleration exerted on the spacecraft which is equivalent to time-dependent 
oscillatory forces that push the bubble in the combined directions of down-and-up 
and leftward-and- rightward as the bubble is rotating with respect to the 
spacecraft rotation axis. Figure 6 exactly illustrate these combined effects of 
equally weighted gravity gradient and jitter accelerations which drive the bubble 
in a manner of combined oscillations of the left-side-up, right- side -down and 
whole bubble down-and-up motions. 

Figure 6 shows the time sequence evolution of the three-dimensional 
behavior of the liquid-vapor interface oscillations driven by equally weighted 
gravity gradient and jitter accelerations. A combined oscillations of left-side 
up and right-side down together with down-and-up bubble oscillations modified by 
leftward- and- rightward oscillations of bubble in the vertical cross-section of 
doughnut profiles indicates very important characteristics of combined gravity 
gradient and jitter accelerations of equal weight acting on the spacecraft when 
it is orbiting around the Earth. 

(III-C) Combined Gravity Gradient and 1CT 6 g 0 Jitter 
Accelerations Acting on Bubble Oscillations 
In this case, combined effects of gravity gradient and jitter accelerations 
are completely dominated by the gravity jitter acceleration when the 10" 6 g c and 
higher gravity jitter acceleration are applied simultaneously with gravity 
gradient acceleration (which is on the order of 10“ 7 g Q as shown in Figure 2) for 
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the excitation of slosh wave along the liquid-vapor interface oscillations. 

Figure 7 shows the time sequence evolution of the dynamical behavior of the 
bubble oscillations driven by gravity jitter-dominated acceleration of 0.1 Hz low 
jitter frequency. It clearly shows that there are a series of asymmetric 
oscillations excited along the surface of liquid-vapor interface driven by 
asymmetric gravity jitter-dominated acceleration. 

Figure 7 shows the time sequence evolution of the three-dimensional 
dynamical behavior of the capillary effect governed interface oscillations driven 
by gravity jitter-dominated acceleration of 0.1 Hz low jitter frequency. It 
clearly shows that there are a series of asymmetric oscillations excited along 
the surface of liquid-vapor interface driven by asymmetric gravity jitter- 
dominated acceleration. 

Careful examination for the case of gravity jitter acceleration- dominate 
liquid-vapor interface oscillations, gravity jitter acceleration, shown in 
Equations (6) and (7) which is also illustrated in Figure 3, indicate that there 
is a sinusoidal oscillation longitudinal to the direction aligned with spacecraft 
mass center to Earth center (parallel to unit vector r c shown in Figure 1) . As 
we indicated in Equation (5), the angle varies with time. This phenomena 
shows that gravity jitter acceleration exerted on the spacecraft is equivalent 
to time -dependent oscillatory forces which push the bubble in the combined 
directions of down-and-up (see z-component of gravity jitter acceleration shown 
in Figure 3) and lef tward-and- rightward (see x-and-y components of gravity jitter 
accelerations shown in Figure 3) as the bubble is rotating with respect to 
spacecraft rotation axis. Down-and-up motion of the bubble is due to the fact 
that z-component gravity jitter acceleration decreases from maximum downward 
direction to minimum downward direction, during the time period between t — 0 to 
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t — 300 s, which push liquid downward and bubble upward; while gravity jitter 
acceleration changes to its downward direction and magnitude from minimum upward 
to maximum upward, and then back to minimum upward, during the time periods 
between t - 300 s to t - 900 s, which push liquid upward and bubble downward; in 
the meanwhile, gravity jitter acceleration changes to its downward direction and 
magnitude during the period from t - 900 s to 1200 s which repeats the half cycle 
story of time period from t - 0 to t - 300 s. Lef tward- and -rightward 
oscillations of bubble are the results of positive and negative directions of 
gravity jitter accelerations in x- and y-components as that illustrated in Figure 
3. Positive direction gravity jitter acceleration implies that liquid is pushed 
rightward and the bubble is driven leftward. Time variations of gravity jitter 
accelerations, shown in Figure 3, characterize the leftward and rightward 
oscillations of bubbles shown in Figure 7. 

In addition to the 0.1 Hz low frequency gravity jitter-dominated 
acceleration, 1.0 Hz medium frequency and 10 Hz high frequency gravity jitter 
acceleration driven liquid-vapor interface oscillations also have been 
investigated. Characteristics of these three ranges of gravity jitter 
frequencies have been well -documented in the earlier studies 3 " 6 . Because of page 
limitation, oscillations driven by the time-dependent variations of the direction 
of gravity jitter-dominated acceleration with medium and high frequencies will 
not be shown in this paper. Results show that lower frequency gravity jitter- 
dominated acceleration contribute more in the driving of asymmetric profiles of 
the interface than that of the higher frequency gravity jitter-dominated 
acceleration. These results are agreeable with the earlier studies*" 7 . 

VI. Mathematical Formulation of the Simulation of Sloshing Dynamics 
Induced Fluid Force and Torque Fluctuations 


13 



To model the simulation of forces and torques exerted on the dewar by large 
amplitude liquid motions, the fluid stresses are decomposed into tangential n t 
and normal components relative to the walls. These expressions are: 




+ 




( 8 ) 




( 9 ) 


where P is the thermodynamic pressure, u a is the fluid velocity in the a 
direction, and t a are unit vectors normal and tangential to the walls, and 5 a 0 
is the Kronecker delta; subscripts a and fi imply the directions. In the 
computation of how fluid force and torque fed back to the container, mathematical 
formulation of non- inertia frame (container bound coordinate) derived earlier 
have to be transformed back to inertia frame (Earth bound coordinate) to carry 
out the calculation 15,28 . Also, to accommodate the pitching, yawing, and rolling 
of the spacecraft, it is convenient here to use Cartesian coordinate (x, y, z) 
with corresponding velocity components. For GP-B, the axis of rotation is always 
fixed at the proof mass located at the mass center of the dewar (x c , y c , z c ) - 
(0, 0, L/2). Detailed derivations and expressions for the fluid-stress forces 
(F x , F y , F z ) and moments (M*, My, M z ) are given in our recent work 15,28 and so will 
not be repeated here. The moment arm relating the forces and moments is given 
by: 
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where L*, Ly, and L z denote the components of the moment arm along x, 
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V. Characteristics of the Simulation of Slosh Waves Induced Fluid 
Force and Torque Fluctuations Exerted on the Dewar Container 
With references to Sections III and IV, one can calculate the fluid force 
and the associated torques exerted on the dewar container of spacecraft 
propellant tank. 

(V-A) Combined Gravity Gradient and 1CT 8 Jitter 
Acceleration Driven Fluid Forces Fluctuations 
As we mentioned earlier, the combined effects of gravity gradient and 
jitter accelerations are completely dominated by gravity gradient acceleration 
when the 10' 8 and less gravity jitter accelerations are applied simultaneously 
with gravity gradient acceleration which drive viscous stress and its moment 
fluctuations . 

Figure 8 shows the computed time variation of the fluctuations of fluid 
forces exerted on the dewar container driven by gravity gradient -dominated 
acceleration. This figure shows the following results: (a) The values of fluid 

force fluctuations are (AF X , AF y , AF Z ) - (23.7, 18.3, 61.7) dynes, it clearly 
indicates AF Z > AF X > AF y . The maximum absolute values of fluid force are Max 
(|F X |, l F y |, |Fj) - (17.4, 11.5, 32.4) dynes. It also indicates |F 2 | > |f x | > 
I F y | . (b) The initial values of F x and F y start from zero value while that of 

the F z starts from non-zero value. This is due to the fact that the shape of 
rotating bubble is symmetric with respect to the x and y axes which makes F x and 
F y equal to zero at time t — 0 while the resultant force along the rotating axis 
of z is not equal to zero because of pre-existing gravity force acting along the 
z-axis at the very beginning, (c) The characteristics of fluid forces and their 
fluctuations are very much following the trend of dominated gravity gradient 
acceleration, shown in Figure 2. In other words, gravity gradient acceleration 
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as a major input force, modulated by sloshing dynamics of the fluid system in the 
dewar, activate on the dewar container with response force. In reality, the 
fluid system, modulated by sloshing dynamics, acts as a damper in the spacecraft 
to damp out the acting force of gravity gradient acceleration. 

Figures 9(A), (B) and (C) show variations of fluid torque exerted on the 
dewar container due to slosh waves driven by gravity gradient-dominated 
acceleration along the x, y and z axes, respectively. The values of fluid torque 
fluctuations are (AM*, AMy, AM Z ) - (761.9, 667.3, 0.06) dyne -cm. The maximum 
absolute values of fluid torque are Max (|M*| f | My | , |M Z |) - (392.9, 343.4, 0.50) 
dyne* cm. It shows AM* > AMy > AM Z and |M*| > | My | > |M Z |. As My, M*, M z are 

determined by the Factors (F x , L z , F z , L*) , (F z , Ly, F y , L.) and (F y , L*, F x , Ly) , 
respectively, it was illustrated in Figures 8 and 10 that the magnitudes of force 
fluctuations give |F Z | > |F X | > |F y | and that of moment arm fluctuations give 
> 1 1* I > M in which time averages of L* and Ly were approaching zero. This 
explains why the moment of M z - L*F y - LyF x - 0 even though there is a rolling 
moment with angular velocity of 0.1 rpm in which the induced rolling moment 
exerted by the viscous rotating fluids is small in comparison with the induced 
pitch and yaw moments due to asymmetric slosh wave excitation. 

Figure 10 shows time fluctuations of moment arms of fluid torques exerted 
on the dewar container due to slosh waves excited by gravity gradient-dominated 
acceleration. The values of moment arm fluctuations are (AL*, ALy, AL Z ) — (22.8, 
19.0, 45.5) cm. The maximum absolute values of moment arm are Max (|L*|, |Ly| 
|L*I> “ (12.2, 10.2, 103.4) cm. It shows AL Z > AL* > ALy and |L Z | > | | > |Ly| . 
This trend is very much similar to the driving force of gravity gradient 
acceleration, shown in Figure 2. 

(V-B) Combined Gravity Gradient and 10~ 7 Jitter 
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Acceleration Driven Fluid Forces Fluctuations 


As we mentioned earlier, the combined effects of accelerations are equally 
weighted by both the gravity gradient and jitter accelerations when 10“ 7 ^ 
gravity jitter acceleration is jointly applied simultaneously with gravity 
gradient acceleration which drive viscous stress and its moment fluctuations. 

Fluctuations of fluid force, fluid torque and moment arm of fluid torques 
driven by gravity jitter acceleration with frequencies of 0.1 , 1.0 and 10 Hz 
have been simulated and investigated. Because of page limitation, results of 
fluctuations driven by 0.1 Hz gravity jitter frequency only will be presented in 
this study. Figure 11 shows the computed time variation of the fluctuations of 
fluid force exerted on the dewar container driven by equally weighted gravity 
gradient and jitter accelerations with a frequency of 0.1 Hz. The values of 
fluid force fluctuations are (AF X , AF y , AF Z ) - (62.4, 57.6, 110.4) dynes and the 
maximum absolute values of fluid force are Max (|F X |, |F y |, | F 2 1 ) - (33.0, 32.3, 
68.1) dynes. It shows AF Z > AF X > AF y and | F 2 1 > |F X | > |F y | . Comparison between 
Figures 8 and 11 for fluid force fluctuations driven by gravity gradient- 
dominated acceleration and equal weight gravity gradient and jitter accelerations 
show that they are in the same order, and it also shows that equal weight gravity 
gradient and jitter accelerations excited values and fluctuations of fluid force 
exerted on the dewar container have a greater value of magnitudes than that 
driven by gravity gradient-dominated acceleration. Comparison among Figures 2,3 
and 11 show that Figure 11 is equally modulated by Figures 2 and 3. In other 
words, the response of sloshing dynamics modulated forces are equally influenced 
by both gravity gradient and jitter accelerations with smaller amplitudes of 
fluctuations due to sloshing dynamics. 

Figures 12(A), (B) and (C) show variations of fluid torques exerted on the 
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container wall due to slosh waves excited by equal weight gravity gradient and 
jitter acceleration with frequency of 0.1 Hz along the x, y and z axes, 
respectively. The values of fluid torque fluctuations are (AM*, AMy, AM Z ) - 
(1038.8, 823.2, 0.07) dyne cm and the maximum absolute values of fluid torque are 
Max (|M*|, iMyl, |M,|) - (538.5, 484.7, 0.50) dyne cm. It shows AM* > AMy > AM* 
and |M*| > |My| > | M z | . Characteristics of the fluctuations of fluid torques 
driven by equal weight gravity gradient and jitter accelerations are similar to 
that drawn for the fluctuations of fluid torques driven by gravity gradient - 
dominated acceleration with greater amplitudes and fluctuations. 

Figure 13 shows time fluctuations of moment arm of fluid torque exerted on 
the dewar container due to slosh waves excited by equal weight gravity gradient 
and jitter accelerations. The values of moment arm fluctuations are (AL*, ALy, 
ALz) - (20.3, 14.8, 33.3) cm. The maximum absolute values of moment arm are Max 
<IM, |Ly|, | l* | ) - (11.5, 9.3, 99.0) cm. It shows AL* > AL* > ALy and |L*| 

> |i,l > |L,I- 

(V-C) Combined Gravity Gradient and 10 _6 g o Jitter 
Accelerations Driven Fluid Forces Fluctuations 

As we mentioned earlier, the combined effects of gravity gradient and 
jitter accelerations are completely dominated by the gravity jitter 
acceleration 16 when the 10’ 6 g Q and higher jitter acceleration are applied 
simultaneously with gravity gradient acceleration which drive fluid force and its 
torque fluctuations . 

Figure 14 shows the computed time variation of the fluctuations of fluid 
force exerted on the dewar container driven by gravity jitter-dominated 
acceleration with a frequency of 0,1 Hz. The values of fluid force fluctuations 
are (AF X , AF y , AF Z ) — (383.0, 315.5, 619.7) dynes and the maximum absolute values 
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of fluid force are Max (|F X |, |F y |, |F Z |) - (192.2, 172.3, 359.5) dynes. It 
shows AF Z > AF X > AF y and | F x | > |F X | > j F y | . Comparison of input and response 
of forces between Figures 3 and 14 show that these two sets of curves are very 
much similar in characteristics and trends with smaller amplitudes and 
fluctuations shown in the response curve. In other words, it' indicates that the 
fluid system, with modulation of sloshing dynamics, acts as a damping machine in 
input and response circuit. Comparison between Figures 8 and 14 for fluid force 
fluctuations driven by gravity gradient -dominated acceleration and gravity 
j itter-dominated acceleration show that gravity j itter-dominated acceleration are 
responsible for the excitation of greater values and magnitudes than that driven 
by gravity gradient-dominated acceleration. 

Figures 15(A), (B) and (C) show variations of fluid torque exerted on the 
container wall due to slosh waves excited by gravity j itter-dominated 
acceleration with frequency of 0.1 Hz along the x, y and z axes, respectively. 
The values of fluid torque fluctuations are (AM*, AMy, AM*) - (1612.2, 1957.9, 
0.10) dyne-cm and the maximum absolute values of fluid torque are Max (|m*|, 
Kl. |Mj) - (1013.7, 1139.3, 0.52) dyne*cm. It shows AMy > AM* > AM Z and |My| 
> K\ > |M,|. Characteristics of the fluctuations of fluid torques driven by 
gravity jitter-dominated acceleration are similar to that drawn for the 
fluctuations of fluid torque driven by gravity jitter-dominated acceleration. 
However, gravity j itter-dominated acceleration excited values and fluctuations 
of fluid torque exerted on the dewar container have a greater value and magnitude 
than that driven by gravity gradient-dominated acceleration. 

Figure 16 shows time fluctuations of moment arm of fluid torque exerted on 
the dewar container due to slosh waves excited by asymmetric gravity jitter- 
dominated acceleration. The values of moment arm fluctuations are (AL*, ALy, 
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AL*) — (4.6 , 3.9, 5.3) cm. The maximum absolute values of moment arm are Max 
(|L*|. |1*|, |L,|) - (2.4, 3.3, 77.8) cm. It shows AL Z > ALy > AL*, and |L Z | > 
IM > IlJ- Comparison of gravity gradient -dominated and gravity jitter- 
dominated accelerations driven fluctuations provide the following conclusions: 
(a) Torsional moraent-and twisting force -equivalent gravity gradient-dominated 
acceleration exerted on the rotating dewar container produces smaller order of 
values and fluctuations of fluid force and fluid torque than that driven by 
gravity jitter-dominated acceleration. (b) Dynamics of bubble (liquid-vapor 
interface) driven by torsional moment-and twisting force -equivalent gravity 
gradient-dominated acceleration produces greater fluctuations and magnitude of 
fluid moment arm than that produced by up and down oscillations of bubble driven 
by gravity jitter-dominated acceleration. (c) Input and response relation 
between combined gravity gradient/jitter accelerations and fluid forces activated 
on the dewar container through the modulation of sloshing dynamics show that the 
characteristics and trends are similar. However, the fluid system acts a damper 
which, in a realistical sense, helps damping out the spacecraft acceleration 
activated by the orbital forces. 

V. Discussion and Conclusion 

Dynamics of spacecraft disturbances activated by the imbalance of liquid 
propellant 29 * 30 driven by gravity gradient-dominated and gravity jitter-dominated 
accelerations, in particular the actuation of the disturbances of spacecraft 
dynamics through the effect of surface tension on partially- filled rotating dewar 
containers (liquid helium and helium vapor) with time-dependent microgravity 
environment have been carried out by numerically computing the three-dimensional 
non- inertia frame Navier-Stokes equations subject to the initial and boundary 
conditions. The initial condition for the liquid-vapor interface profiles were 
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adopted from the steady- state formulation of rotating dewar developed by Hung and 
Leslie 17 . 

In this study, spacecraft gravity turn-around time is assumed to be 1200s. 
The equivalent twisting force and torsional moment of gravity gradient -dominated 
acceleration, acted on the rotating dewar orbiting around the Earth, induce left- 
side-up and right -side -down in the horizontal cross section of doughnut profiles. 
For the case of bubble oscillations driven by gravity j itter-dominated 
acceleration, up-and-down and leftward- and- rightward motion have been resulted. 

Fluid force and fluid torque moment exerted on the dewar wall along the 
yawing and pitching axes caused by slosh wave excitation driven by gravity 
gradient-dominated and gravity jitter-dominated accelerations in partially 
liquid- filled container, in addition to the existing rotating motion along the 
rolling axes, have been simulated and investigated. Results show that there are 
a series of large amplitude fluctuation of fluid forces, fluid torque and moment 
arms exerted on the dewar container of spacecraft propellant tank which in turn 
disturb the dynamics of spacecraft along yawing and pitching axes. In 
particular, slosh waves excited by the asymmetry gravity gradient-dominated and 
gravity jitter-dominated accelerations contribute greatly to the large amplitude 
fluctuations of fluid forces and fluid torque exerted on the container and the 
dynamics of spacecraft along the pitching and yawing axes, in addition to the 
existing rolling axis of spacecraft 29,30 . 

Time -dependent variations in the direction of gravity gradient-dominated 
and gravity jitter-dominated accelerations imposed on the spacecraft will change 
and modify the asymmetry distribution of liquid-vapor in the rotating dewar 
container and the dynamic system of spacecraft 31 . This means that the dynamic 
system of spacecraft through the fluid forces and fluid torque exerted on the 
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container wall will be deeply disturbed by the presence of asymmetric 
fluctuations in liquid-vapor interface 29 ' 30 . In this study, the differences in 
gravity gradient-dominated acceleration and gravity j itter-dominated acceleration 
induced slosh waves and how these slosh wave excitation affected the dynamic 
system of spacecraft through the fluctuations in fluid force and fluid torque 
exerted on the container wall distribution have been simulated and studied. 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Coordinate system for the computation of gravity gradient 
acceleration, 0 g is the location of geometric center of cylinder. 
Direction of r c points toward Earth center. (x' , y' , z’) denotes 
absolute coordinate system, while (x, y, z) illustrates non-inertia 
frame spacecraft bound coordinate system. 

Time variation of gravity gradient acceleration acting on fluid mass 
located at (r, 9, z) - (40 cm, n/b, 10 cm) for turn-around period of 
1200 s with rotating speed of 0.1 rpm. (A) Along x-direction; (B) 
Along y-direction; (C) Along z-direction. 

Time variation of gravity jitter acceleration acting on the fluid 
system with a value of 10‘ 8 , 10’ 7 and 10‘ 6 g Q , rotating speed of 0.1 
rpm, turn-around period of 1200 s and gravity jitter frequency of 0.1 
Hz. (A) Along x-direction; (B) Along y-direction; (C) Along z- 
direction. 

Initial profiles of liquid-vapor interface for Gravity Probe-B module 
of rotating dewar with gravity of 10‘ 7 g 0 , rotating speed of 0.1 rpm 
and direction of background gravity at - 0°. (A) In r-z plane at 

* “ 0° and 180°, (B) In r-0 plane at 6 - 90° and 270°, (C) In r -9 

plane at z - 108 cm, and (D) Three-dimensional liquid-vapor interface 
profile . 

Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by combined gravity gradient 
and jitter accelerations with gravity of 10" 8 g Q . Rotating speed of 
dewar is 0.1 rpm, and the period of gravity direction turn-around 
time is 1200 s. 
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Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 10 


Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by combined gravity gradient 
and jitter accelerations of 0.1 Hz frequency, gravity of 10‘ 7 ^ with 
the period of gravity direction turn-around time - 1200 s, and 
rotating speed of 0.1 rpm. 

Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by combined gravity gradient 
and jitter accelerations of 0.1 Hz frequency, gravity of 10~ 6 g^, with 
the period of gravity direction turn-around time - 1200 s, and 
rotating speed of 0.1 rpm. 

Time sequence of the fluctuations of fluid stress forces exerted on 
the dewar container due to the slosh wave excitation driven by 
combined gravity gradient and jitter accelerations with 10' 8 ^ under 
rotating speed of 0.1 rpm and gravity turn-around time of 1200 s. 
(A) Along x-direction; (B) Along y-direction; and (C) Along z- 
direction. 

Time sequence of the fluctuations of fluid stress moment exerted on 
the dewar container due to the slosh wave excitation driven by 
combined gravity gradient and jitter accelerations with 1O -0 g 0 under 
rotating speed of 0.1 rpm and gravity turn-around time of 1200 s. 
(A) Along x-direction; (B) Along y-direction; and (C) Along z- 
direction. 

Time sequence of the fluctuations of moment arm of viscous 
stress moment exerted on the dewar container due to the slosh wave, 
excitation driven by combined gravity gradient and jittur 
accelerations with 10' 8 g Q under rotating speed of 0.1 rpm and 
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Figure 11 


Figure 12 


Figure 13 


Figure 14 


gravity turn-around time of 1200 s. (A) Along x-direction (pitching 
moment); (B) Along y-direction (yawing moment) ;and (C) Along z- 
direction (rolling moment). 

Time sequence of the fluctuations of fluid stress moments exerted on 
the dewar container due to the slosh wave excitation driven by 
combined gravity gradient and jitter accelerations with 10' 7 g^ 
rotating speed of 0.1 rpm, gravity turn-around time of 1200 s and 
gravity jitter frequency of 0.1 Hz. (A) Along x-direction (pitching 
moment) ; (B) Along y-direction (yawing moment) ; and (C) Along z- 
direction (rolling moment) . 

Time sequence of the fluctuations of fluid stress moments exerted on 
the dewar container due to the slosh wave excitation driven by 
combined gravity gradient and jitter accelerations with 10" 7 g 0 
rotating speed of 0.1 rpm, gravity turn-around time of 1200 s and 
gravity jitter frequency of 0.1 Hz. (A) Along x-direction fluid 
stress moment; (B) Along y-direction fluid stress moment; and (C) 
Along z-direction fluid stress moment. 

Time sequence of the fluctuations of moment arm of fluid stress 
moments exerted on the dewar container driven by combined gravity 
gradient and jitter accelerations with 10' 7 g^ rotating speed of 0.1 
r P m » gravity turn-around time of 1200 s and gravity jitter frequency 
of 0.1 Hz. (A) x-direction moment arm; (B) y-direction moment arm; 
and (C) z-direction moment arm. 

Time sequence of the fluctuations of fluid stress moments exerted on 
the dewar container due to the slosh wave excitation driven by 
combined gravity gradient and jitter accelerations with 10‘ 7 g Q , 
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rotating speed of 0.1 rpm, gravity turn-around time of 1200 s and 
gravity jitter frequency of 0.1 Hz. (A) Along x-direction (pitching 


moment); (B) Along y-direction (yawing moment); and (C) Along z- 
direction (rolling moment). 

Figure 15 Time sequence of the fluctuations of fluid stress moments exerted on 
the dewar container due to the slosh wave excitation driven by 
combined gravity gradient and jitter accelerations with 10" 6 g ot 
rotating speed of 0.1 rpm, gravity turn-around time of 1200 s and 
gravity jitter frequency of 0.1 Hz. (A) Along x-direction fluid 
stress moment; (B) Along y-direction fluid stress moment; and (C) 
Along z -direction fluid stress moment. 

Figure 16 Time sequence of the fluctuations of moment arm of fluid stress 
moments exerted on the dewar container driven by combined gravity 
gradient and jitter accelerations with 10 -6 g 0f rotating speed of 0.1 
r P m > gravity turn-around time of 1200 s and gravity jitter frequency 
of 0.1 Hz. (A) x-direction moment arm; (B) y-direction moment arm; 
and (C) z- direction moment arm. 
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INITIAL PROFILES OF LIQUID HELIUM AND VAPOR INTERrACE 
Rotating Dewar Under the Effect of Gravity Jitter 
g = 10" 7 gQ,^= 0.1 rpm , / Vi« 0 , f = 0.1 Hz 
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(D) Three-Dimensional Profiles 
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